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Abstract. For gradiometric three-Josephson-j unction flux qubits, we perform a 
systematic study on the tuning of the minimal transition frequency, the so-called qubit 
gap. By replacing one of the qubit's Josephson junctions by a dc SQUID, the critical 
current of this SQUID and, in turn, the qubit gap can be tuned in situ by a control flux 
threading the SQUID loop. We present spectroscopic measurements demonstrating a 
well-deflned controllability of the qubit gap between zero and more than 10 GHz. In 
the future, this enables one to tune the qubit into and out of resonance with other 
superconducting quantum circuits, while operating the qubit at its symmetry point 
with optimal dephasing properties. The experimental data agree very well with model 
calculations based on the full qubit Hamiltonian. From a numerical fit, we determine 
the Josephson coupling and the charging energies of the qubit junctions. The derived 
values agree well with those measured for other junctions fabricated on the same chip. 
We also demonstrate the biasing of gradiometric flux qubits near the symmetry point 
by trapping an odd number of flux quanta in the gradiometer loop. In this way, 
we study the effect of the significant kinetic inductance, thereby obtaining valuable 
information for the qubit design. 



PACS numbers: 03.67.Lx, 74.50.+r, 85.25.Cp 
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1. Introduction 

Superconducting quantum circuits are promising for the implementation of solid 
state quantum information systems [l|[2] and the realization of fascinating quantum- 
optical experiments in the microwave regime (3||4]. In particular, the coupling of 
superconducting quantum two-level systems (qubits) to microwave resonators has been 
successful |5-7|, resulting in the rapid development of the prospering field called 
circuit quantum electrodynamics (QED). In circuit QED, strong [5||6] and ultrastrong 
coupling [8- lOj of superconducting quantum bits to the electromagnetic modes of high 



quality factor microwave resonators has been demonstrated. Circuit QED also has been 



used to generate non-classical states of light 11, 12 , to establish single artificial atom 



masing |13|, to realize controlled symmetry breaking (14,15 , or to implement quantum 
gates and algorithms n!6,17 



Nowadays, the most popular superconducting qubits are the cooper pair box (5|[18| 



the transmon qubit 12,16,19 , the phase qubit 11 20 21 , the fiuxonium 22, and 



the persistent current or fiux qubit j^, 23-25 . For the implementation of circuit QED 
experiments, transmon and phase qubits have been used most often for several reasons. 
First, the relevant qubit parameters can be controlled within sufl&ciently narrow margins 
in the fabrication process. Second, a controlled coupling/ decoupling to a microwave 
resonator acting as a quantum bus is possible by a fast change of the qubit 's transition 
frequency. Third, the coherence properties of the qubit do not strongly degrade during 
such operations. Unfortunately, the original design of the fiux qubit f23| consisting 
of a superconducting loop intersected by three Josephson junctions cannot fulfill these 
requirements simultaneously. First, although the fiux qubit 's transition frequency cjq 
can be varied over a wide range by applying an external magnetic fiux, the coherence 
time of the fiux qubit rapidly decreases when tuning the qubit away from its symmetry 
point with minimum transition frequency cjq = A. Only at its symmetry point the 
fiux qubit is well protected from the relevant l//-noise and coherence times exceeding 
1 /is can be reached fl]. Second, for the fiux qubit the minimal energy splitting hA 
between the ground and excited state depends exponentially on the critical current 
Ic and capacitance Cj of the Josephson junctions |25| and therefore is difiicult to 
precisely control in fabrication. This does not allow to fabricate fiux qubits with well- 
defined A values, which are, for example, close to the resonant angular frequency 
of superconducting microwave resonators. On the other hand, fiux qubits have specific 
advantages. First, the anharmonicity of fiux qubits, that is, the separation of the excited 
state from the third level relative to cjq, is by one or two orders of magnitude larger than 
for transmon and phase qubits, allowing for fast qubit operations without leakage to 
higher states. Second, fiux qubits can be coupled ultrastrongly to resonators. Relative 
coupling strengths g/cjUj; > 0.1, where is the resonator frequency and g the coupling 
strength, have been demonstrated |8,9|. The reason is that the coupling is inductive 
for the fiux qubit and capacitive for the transmon/phase qubit, respectively. As the 
inductive and capacitive coupling is proportional to $o^rms and 2el^nis5 respectively. 
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the ratio of the couphng strengths is (/i/4e^)(/rms/Kms) = (/i/4e^)/Zr ^ 10. Here, 
$0 = ^/2e is the flux quantum, e the electron charge, /rms and Kms the zero-point 
current and voltage fluctuations in the resonator and the characteristic impedance 
of the resonator. 

In order to overcome the drawback of flxed minimum energy splitting /lA in 
superconducting flux qubits, Orlando et al [25] proposed a modifled flux qubit 
design, which subsequently has been implemented by Paauw et al |26| and meanwhile 



successfully used in several experiments flO, 27-30 . In this tunable-gap flux qubit 



one of the Josephson junctions, the so-called a-junction, is replaced by a small loop 
with two Josephson junctions. This dc superconducting quantum interference device 
(SQUID) acts as a junction whose critical current can be controlled by the flux threading 
the SQUID loop. As a consequence, also the qubit gap A can be tuned in such a 
conflguration. The additional control allows for a fast variation of the qubit transition 
frequency cjq, while operating the flux qubit at its symmetry point where the coherence 



properties are optimum [10,26,28 . In this article, we report on the fabrication and 
systematic study of tunable-gap gradiometric flux qubits, following the design proposed 
in We show that the energy splitting at the symmetry point can be varied in a 

controlled way between values close to zero and A/27r > 10 GHz. In particular, our 
analysis extends to the case of multiple flux quanta trapped in the gradiometer loop. 
In this way, we obtain detailed insight into the tuning mechanism. In section [2| we 
first introduce the foundations of fixed-gap and tunable-gap fiux qubits required for the 
analysis of the experimental data. In particular, we discuss the possibility of fiux-biasing 
gradiometric fiux qubits at the symmetry point by freezing in an odd number of fiux 
quanta during cool-down and the effect of the significant kinetic inductance of the narrow 
superconducting lines of the qubit loop. After briefiy introducing the experimental 
techniques in section [3| we present the experimental data and their analysis in section [4] 
before concluding with a brief summary in section [5| 



2. The flux qubit 

In the following we briefiy summarize the foundations of fixed-gap and tunable-gap 
persistent current or fiux qubits as well as their gradiometric versions. We derive the 
relevant expressions used in the evaluation of our experimental data. 



2.1. The fixed-gap flux qubit 

The simplest version of the fiux qubit [cf. figure[l](a)] consists of a small superconducting 
loop with a diameter of the order of 10 /xm intersected by three Josephson junctions (J J) 
with lateral dimensions of the order of 100 nm [23]. While two of these J J have the same 
area (typically, Aj 0.03 /xm^ in our experiments) and, hence, the same critical current 
(typically, Ic — 600 nA), the third J J, the so-called a-junction, has a reduced area 
Aa = aAj with a 0.6 — 0.8, resulting in a reduced critical current Ic^a = and 
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reduced junction capacitance Ca = cxCj. Since a = aois fixed in the fabrication process, 
the qubit gap A is also fixed. Consequently, this version of the flux qubit is called fixed- 
gap flux qubit. For a ^ 0.6 — 0.8, the two-dimensional potential energy landscape 
of the flux qubit can be simplified. At the symmetry point, where the magnetic fiux 
through the loop is equal to (n + |)$05 with n being an integer, the potential can be 
reduced to a one-dimensional double well |25|. The two minima of this potential are 
associated with two degenerate persistent current states, corresponding to clockwise and 
counter-clockwise circulating persistent currents ±/p. Due to the finite tunnel coupling 
of these states their degeneracy is lifted. The resulting symmetric and anti-symmetric 
superposition states form the ground and excited state of the fiux qubit separated by 
the minimal energy splitting HA. In the basis of the persistent current states and near 
the symmetry point, the Hamiltonian describing the fiux qubit can be written as |25j 

n = hea, - hAa, . (1) 

Here, and ax are the Pauli spin operators, he = 2/p5$ is the magnetic energy bias, 
and 5$ = $o[/ ~ {ti + |)] the deviation of the fiux $ threading the loop from a half- 
integer multiple of $o- The quantity / = $/$o is the magnetic frustration of the qubit 
loop and n an integer. The transition frequency between ground and excited state can 
be written as 

CJq = VS^ + A2 (2) 

and the qubit gap A becomes [25j 



HEAHa^-D /4„(i+2a)t| . (3) 



a{l + 2a) I ' ' V E, 

Here, Ej = hl^/2e is the Josephson coupling energy, = /2C the charging energy 
and a{a) = ^1 — {l/Aa^) — [arccos(l/2a)/2a] with a{a) 0.15 for a = 0.7. Obviously, 
A is determined by the critical current /c and the capacitance Cj of the Josephson 
junctions as well as by a = a^. All these parameters are fixed by the fabrication 
process. Decreasing a from 1 to 0.5 results in a strong increase of the exponential 
factor. At the same time, the prefactor (attempt frequency) decreases from the plasma 
frequency of the JJ to zero, because the double well potential becomes a single well at 
a = 0.5. Since the exponential factor dominates within the major part of the interval 
0.5 < a < 1, a strong increase of A is obtained by reducing a. 

It is evident from ([T]) and ^ that it is possible to tune cjq by varying either e 
or A. Varying the energy bias s is simply achieved by changing via the magnetic 
field generated by an external solenoid or a current fed through an on-chip control 
line. However, 5^ causes a shift of the qubit operation point away from the anti- 
crossing point with minimal transition frequency A. As the energy of the fiux qubit is 
stationary with respect to small variations of the applied magnetic fiux {dUci/dS^ = 0) 
only for / = + |), any shift away from this symmetry point makes the fiux qubit 
more susceptible to magnetic fiux noise and significantly deteriorates the coherence 
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Figure 1. Circuit schematics of (a) the three- Josephson-junction (3- J J) flux qubit 
with a-j unction, (b) the simplest 3- J J flux qubit with the tunable a-j unction realized 
by a dc SQUID, (c) the gradiometric 3-JJ flux qubit without tunable a-junction, and 
(d) the gradiometric 3-JJ flux qubit with tunable a-j unction. The a- and e-lines can 
be used to change the magnetic frustration of the a- and qubit loop independently. 



properties |3Tp2] . Since the fast tuning of cjq of flux qubits is a prerequisite for numerous 
circuit QED experiments [33|[34j , it is desirable to realize a tuning of cjq by a variation 
of A. 

2.2. The tunable- gap flux qubit 

According to ([2]), cjq cannot only be tuned by varying s but also by varying A. This is 
advantageous, since the operation point of the qubit stays at the symmetry point with 
optimal coherence properties. Flux qubits with tunable A are called tunable- gap flux 
qubits. As pointed out by Paauw et al. (26|, an in situ tunability of A is achieved by 
replacing the a-junction by a small a-loop containing two J J [cf. figure [l](b)]. That is, 
the a-junction is replaced by a dc SQUID. Then, the critical current Ic^a of the a-loop 
and, in turn, the qubit gap A can be tuned by a control flux threading the a-loop. 
If we choose the area of the junction in the a-loop to O.baoAj^ we obtain Ic^a = with 
a = ao| cos(7r$c^/$o)|- Successful implementations of this design have been reported 
recently [35|]38|. 

The magnitude of /p can be calculated to |25| 



with A{a) according to ([3]). We see that for a 0.5 the persistent current /p approaches 
zero. 




(4) 



yielding the a-dependent transition frequency 




(5) 
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2.3. The gradiometric flux qubit 

While the replacement of the a-j unction by an a-loop allows for a tunable qubit gap A, 
applying any flux to the a-loop at the same time changes the flux threading the qubit 
loop and hence the energy bias s of the flux qubit. This is unintenional and has to be 
compensated. To keep the energy bias of the flux qubit constant during variations of 
$CK, a gradiometric design can be used. The gradiometric versions of a flxed-gap and 
tunable-gap flux qubit are shown in flgures[l](c) and (d), respectively. 

For the tunable-gap gradiometric flux qubit of flgure[l](d), an applied homogeneous 
magnetic fleld changes $c^, but does not affect the energy bias of the flux qubit, 
since the screening currents in the two subloops of the eight-shaped gradiometric 
loop cancel each other on the central line. The immediate consequence is that an 
inhomogeneous magnetic fleld is required to adjust the energy bias s of the flux qubit. 
This inhomogeneous fleld can be generated by feeding a small current through the 
so-called 6:-flux line, which couples asymmetrically to the qubit loop [cf. flgure [l](d)]. 
Furthermore, the outer ring of the gradiometric qubit, denoted as the trapping loop, 
can be used to trap an integer number of magnetic flux quanta, e.g. by cooling down 
below Tc in an applied magnetic fleld. This allows for a pre-biasing of the qubit near the 
symmetry point. We note, however, that the exact amount of flux threading the qubit 
loop and the a-loop, respectively, depends on the ratio of the kinetic and geometric 
inductances. Since an understanding of this point is important for a controlled design 



of a gradiometric qubit with tunable gap, it is discussed in more detail in section |2.5 
Within this work, we investigate both flxed-gap and tunable-gap gradiometric 3-JJ flux 
qubits. The former is an ideal model system to study the principle of flux biasing. 



2.4' Flux biasing 

In this subsection we briefly address the flux biasing of gradiometric flux qubits by the 
trapping of magnetic flux in its outer loop, the so-called trapping loop [cf. flgures [l](c) 
and (d)]. Flux biasing is based on the phase coherence of the superconducting state. 
The phase 9 of the macroscopic wave function describing the superconducting state is 
allowed to change only by integer multiples of 27r along a closed integration path: 

V9-ds = 27m . (6) 

In multiply connected superconductors this leads to the expression for the fluxoid 
quantization: 

^ fioXlJs 'ds + J B ' dF = n^o . (7) 

Here, Al is the London penetration depth, /xq the vacuum permeability, F a closed 
integration path encircling the area F, Jg the supercurrent density along F, and B 
the magnetic flux density. The second term on the left hand side represents the total 
magnetic flux $ threading the area F. For superconductors with cross-sectional area 
large compared to the London penetration depth Al, the flrst term vanishes, since 
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one always can find an integration path deep inside the superconductor where the 
supercurrent density Jg = 0. This leads to the expression for flux quantization 

^ = /tr = n , (8) 

saying that the total magnetic flux in a closed superconducting loop such as the trapping 
loop is quantized in units of $o- 

The phase of the superconducting order parameter changes by 27rn around the 
closed trapping loop. Therefore, in the fully symmetric gradiometric qubit designs of 
figures [l](c) and (d) the trapping of an odd number (2n + 1) of flux quanta in the 
trapping loop leads to a phase difference of (2n + l)7r between the points A and B. This 
corresponds to a flux bias of (n + |)$0 5 i-e. a flux bias at the symmetry point. The 
biasing with trapped flux has the advantage that it is not affected by the finite noise of 
current sources required for the biasing with an external magnetic field. On the other 
hand, once a specific flux state has been frozen in, it can no longer be changed without 
heating the sample above Tc. Therefore, in practice flux trapping is often used for pre- 
biasing at an operation point, while an additional magnetic field is used for making fast 
changes around this operation point. In order to enable such flux control the width of 
the superconducting line forming the trapping loop has to be made small enough (of 
the order of Al) to allow for partial penetration of the applied magnetic field. In this 
case, the first term on the left hand side of ([7| becomes relevant. This term is related 
to the kinetic energy of the superconducting condensate or, equivalent ly, the kinetic 
inductance Lk, whereas the second term is related to the field energy or, equivalent ly, 
the geometric inductance Lg of the trapping loop. 

2. 5. The effect of the kinetic inductance 

We next discuss the influence of the kinetic inductance L^, which is no longer negligible 
compared to the geometric inductance Lg of the trapping loop when the width of the 
superconducting lines is reduced to values of the order of Al- In this case the first term 
on the left hand side of ([T]) is no longer negligible. With the supercurrent Idr = JgS 
circulating in the trapping loop, we can rewrite this term as 

lioXlJs ■ds = ^ I,,,e = L^hs, . (9) 

Here, we have introduced the kinetic inductance = /iioXKi/S) of the trapping loop, 
with i its circumference and S its cross-sectional area. With $k = ^k^cir and splitting 
up the total flux $ into a part $ex due to an external applied field and a part $g = L^Idr 
caused by the circulating current in the trapping loop with geometric inductance Lg, 
the fluxoid quantization condition ([7| reads as 
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Introducing the parameter (5 = L^/L]^^ we obtain the expression for the net magnetic 
frustration of the trapping loop to 

1 ^ P 



tr.net 



/ex + fg 



/ex + 



n 



(11) 



The net magnetic frustration of the a-loop in first approximation is obtained by 
multiplying with the area ratio Aa/At^ of the a- and the trapping loop: 

Ar 



fa. 



net 



A Jtr,net 



(12) 



Here, we have neglected effects arising from the fact that the a-loop is not centered in 
the trapping loop. If the geometric inductance is negligible (/3 <C 1), the contribution 
of the circulating screening current is negligible and /tr,net — /ex- In this case, the 
superconducting lines cannot screen magnetic fields and we can change the magnetic 
frustration of the a- and the trapping loop continuously by varying the applied magnetic 
field. In contrast, if the geometric inductance is dominant (/3 1), the screening is so 
strong that we can no longer change the fiux in the loop by varying the applied field. 
The frustration of the trapping loop is fixed to the value /tr,net — frozen in during 
cool-down. This means that also the frustration of the a-loop can no longer be changed 
continuously as desired. 

With the net magnetic frustration (12), the critical current of the a-loop is obtained 

to /c,a(/a,net) = Q^(/a,net)^c with 



<^(/a,net) = C^O I COS (vT^^net)! 



cos TT 



A 



tr 



1 



/ex + 



f3 



n 



(13) 



_l + /5'^'" 1 + /5 

For a suitable width of the superconducting lines and, hence, a suitable value of /3, we 
can vary a both by changing /ex via an external magnetic field and by changing the 
number n of fiux quanta frozen into the trapping loop during cool-down. For example, 
n could be used for pre-biasing at a specific a value and the external magnetic field 
provided by a current sent through an on-chip control line for small variations around 
this value. The pre-biasing with trapped fiux has the advantage that it is not aflFected 
by the noise added by the current source, while the variations with the on-chip control 
line can be very fast. 



For zero applied magnetic field, (13) reduces to 

A. 



<^(/a,net)|/ex=0 = O^Q 



COS TT 



n 



(14) 



Ar 1+13 

This expression applies to the experimental situation, where an odd number (2n + 1) 
of fiux quanta is frozen into the trapping loop to bias the gradiometric fiux qubit at its 
symmetry point and no additional external magnetic field is applied. Fixing — 1 by 
the fabrication process, we can change the number of trapped fiux quanta to choose a 
in the desired regime 0.5 < a < 1. Of course, fiux trapping only allows for a step- wise 
variation of a. For continuous and fast variations of a, magnetic fields generated by 
external solenoids or on-chip control lines have to be used. For a typical value of /3 0.8^ 
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we obtain fa,net — 0.08 n for /ex = 0. This shows that we need only a smaU number of 
trapped flux quanta to significantly modify a. Furthermore, we obtain /tr,net = 0.55/ex 
for n = 0, meaning that about half of the applied magnetic fiux is shielded by the 
trapping loop. 



2.6. The gradiometer quality 

A perfect gradiometer should be completely insensitive to a homogeneous magnetic field. 
However, in reality there are always imperfections such as slight differences of the areas 
Ai and A2 of the two subloops of the eight-shaped gradiometer and/or of the geometric 
inductances Lgi and Lg2 and kinetic inductances Lki and Lk2 of the superconducting 
lines forming the subloops. Due to these imperfections there will be a finite imbalance 



^fimh of the magnetic frustration of the two subloops. According to (10), 5/imb can be 
expressed as 

y^o ^0 ^0 

~ $0 ^ $0 $0 ■ 

With /cir = (n — /ex)$o/(-^g + L]^) we can rewrite this expression to 

<^/imb = /ex r , t + " / , . ■ (16) 



A L„ + L], J \ L„ + L 



— 1/Qgrad,ex — l/Qgrad,n 

The total gradiometer quality Q is given by = (5grid,ex + Qgmd,n- The first term 
describes imbalances of the frustration when a homogeneous external field is applied, 
and the second those when an integer number of fiux quanta is frozen in. Obviously, 
the higher the Q the lower is (^/imb- As shown below, Q values of the order of 500 are 
feasible. 

We note that the ^-fiux line shown in figure [l](c) generates different fiux densities 
B in the two subloops of area yltr,! — Ar,2, leading to different amounts of total fiux 
$1 = r. BdA and $9 = [a BdA. This results in the magnetic frustration 

by the 6:-fiux line, which is used to change the energy bias £ of the gradiometric fiux 
qubit. Correspondingly, the deviation of /12 from the value (n + |) at the symmetry 
point is 

5/12 = /i2- (n + l) . (18) 



2 

3. Experimental techniques 



The fiux qubits used in our study are based on Al thin film structures and Al/AlO^^/Al 
Josephson junctions fabricated by electron beam lithography and two- angle shadow 
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Figure 2. Optical micrograph of the sample chip containing a tunable-gap 
gradiometric flux qubit with the biasing lines and filter structures. The large inset 
shows an enlarged view of the region marked with the blue rectangle: the gradiometric 
qubit with readout dc SQUID as well as the a- and e-flux lines used for tuning the 
frustration of the a- and the qubit loop. The two small insets show scanning electron 
microscopy (SEM) images of a regular (top) and an a-junction (bottom). The position 
of these junction are marked with yellow rectangles in the large inset. 



evaporation on thermally oxidized silicon wafers. Details of the fabrication process can 
be found in [39] . Figure [2] shows optical and SEM micrographs of a sample chip with a 
tunable-gap gradiometric flux qubit. The qubit is surrounded by the readout dc SQUID 
and the control circuitry for the energy bias (6:-flux line) and the qubit gap A (a-flux 
line). The insets show an enlarged view of the qubit region as well as SEM images 
of regular junctions and one of the junctions of the a-loop, which has a reduced area 
0.5ao^j with ao 1. The area of the regular junctions is Aj 0.03 /xm^. 

All measurements have been performed in a dilution refrigerator with a base 
temperature of 30 mK. The qubit state is read out via a dc SQUID inductively coupled 
to the flux qubit (3T,39 . Qubit transitions between the ground and excited state can be 



induced by microwave radiation supplied via an off-chip antenna. The continuous-wave 
(CW) microwave signal applied in our experiments is strong enough to saturate the 
qubit, leading to a 50% population of ground and excited state. The trapping of flux 
quanta in the trapping loop of the gradiometric qubit is obtained by cooling down the 
circuit into the superconducting state in the presence of an appropriate magnetic field. 
Moreover, the sample can be heated up above Tc by applying a suitable heating current 
to an external heater located near the sample. 
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4. Experimental results and discussion 
^.7. Basic parameters 

The critical current density of the Josephson junctions is determined to Jc(30mK) = 
1.5 — 3.5kA/cm^ by measuring the current-voltage characteristics (IVCs) of the readout 
SQUIDs fabricated on the same chip and by determining the junction area Aj by 
scanning electron microscopy. The Jc values can be varied by changing the oxidation 
process. For = 2kA/cm^ and the typical junction area of Aj 0.03 /xm^ we have 
/c 600 nA. The specific capacitance of the junctions is derived from the analysis of 



resonances in the IVCs of the readout SQUIDs 40 . For junctions with = 2kA/cm^, 
we find C/A = (195 ± 10)fF//xm2, resulting in Cj 6fF for Aj ^ 0.03 /xm^. The 
geometric inductance Lg of the superconducting loops are estimated according to |41|. 
In order to estimate the kinetic inductance of the superconducting lines we use the dirty 



limit expression Lk = hp^i/nAoS [42,43 , where Aq = O.lSmeV is the zero temperature 
energy gap of Al. The use of this expression is justified, since the mean free path in 
our 90 nm thick Al films is limited by the film thickness and therefore is much smaller 
than the coherence length ^ 1.5 /xm of Al. The normal resistivity pn is determined 
by suitable test structures fabricated on the same chip. For the cross-sectional area 
S = 500 X 90 nm^ of the superconducting line forming the trapping loop we obtain a 
kinetic inductance per unit length of L]^/i ^ lpH//xm. 

4^2. The fixed- gap flux qubit 

We first discuss the properties of fixed-gap, non-gradiometric fiux qubits serving as 
reference samples. The qubit gap A and the persistent current /p are determined by 



qubit spectroscopy 31,39]. Figure [3| shows typical spectra obtained for two 3- J J fiux 
qubits with fixed a-j unction by sweeping the qubit frustration 5f = / — (^+|) = S^/^q 
at fixed microwave frequency. The qubit state is read out repeatedly by the readout 
dc SQUID. Only at those Sf values where the microwave driving is resonant with the 
qubit transition frequency cjq, a 50% population of the excited state is detected. This 
manifests itself in characteristic peak and dip structures in the switching current /gw of 
the readout SQUID at frequency dependent 5f values. Plotting these values versus the 
microwave frequency as shown in figure [3] yields (jJc^{Sf). Assuming that Jc has the same 
value for all three junctions, the value oi a = = Ac^/Aj can be determined from the 
measured area ratio. Then a two-parameter fit of ([s]) to the spectroscopy data yields A 
and /p = he/ 25^. The spectra in figure [3] are obtained for two fiux qubits diflFering only 
in their values. For = 0.75 and 0.6, we obtain A/27r = 1.39 GHz and 10.76 GHz 
and /p = 583 nA and 283 nA, respectively. Obviously, for values closer to 0.5 (1.0) 
large (small) A and small (large) /p values are obtained in agreement with ^ and A 
consistency check can be made by calculating the /p values from Q. Here, the unknown 
critical current 1^ = JqA is estimated from the measured junction area and using the 
Jc value of the junctions of the readout SQUID. We obtain /p = 619 nA and 310 nA in 




Figure 3. Transition frequency (jJc^/27i plotted versus 8f = / — (n + ^) = S^/^o 
for two fixed-gap flux qubits with (a) ao = 0.75 and (b) ao = 0.6. Also shown is a 
two-parameter fit of the data (black lines) yielding A/27r and /p and the result of a 
numerical simulation based on the diagonalization of the full qubit Hamiltonian. In 
(b) the result of the two-parameter fit and the simulation are almost indistinguishable. 



good agreement with the values derived from the spectroscopy data.rence samples. The 



qubit gap A and the persistent current /p are determined by qubit spectroscopy 31,39 



Figure [3] shows typical spectra obtained for two 3- J J flux qubits with flxed a-junction by 
sweeping the qubit frustration 5/ = / — (n + |) = 5$/$oat flxed microwave frequency. 
The qubit state is read out repeatedly by the readout dc SQUID. Only at those 5f 
values where the microwave driving is resonant with the qubit transition frequency cjq, 
a 50% population of the excited state is detected. This manifests itself in characteristic 
peak and dip structures in the switching current /gw of the readout SQUID at frequency 
dependent 6f values. Plotting these values versus the microwave frequency as shown 
in flgure |3] yields uJc^{Sf). Assuming that Jc has the same value for all three junctions, 
the value of a = ao = A^/Aj can be determined from the measured area ratio. Then 
a two-parameter flt of ^ to the spectroscopy data yields A and /p = he/ 25^. The 
spectra in flgure |3] are obtained for two flux qubits diflFering only in their values. For 
Q^o = 0.75 and 0.6, we obtain A/27r = 1.39 GHz and 10.76 GHz and /p = 583 nA and 
283 nA, respectively. Obviously, for values closer to 0.5 (1.0) large (small) A and 
small (large) /p values are obtained in agreement with ^ and ^ . A consistency check 
can be made by calculating the /p values from Q. Here, the unknown critical current 
/c = JcA is estimated from the measured junction area and using the value of the 
junctions of the readout SQUID. We obtain /p = 619 nA and 310 nA in good agreement 
with the values derived from the spectroscopy data. 

We have also performed numerical simulations based on the diagonalization of the 
full qubit Hamiltonian using and a = ^c^/^j as input parameters. They are based 
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Figure 4. (a) Circuit schematics of the fixed-gap gradiometric flux qubit with readout 
dc SQUID and e-flux hue. The outer loop of the flux qubit (broken ohve fine) forms 
the trapping loop, (b) Scanning electron microscope image of the implemented circuit. 
The inset shows an enlarged view of the a-junction. (c) Transition frequency co'q/27r 
plotted versus Sf = f — {n-\-^) = 6^/^o for a flxed-gap gradiometric flux qubit with 
ao = 0.65. Also shown is the result of a numerical simulation (red line) based on the 
diagonalization of the full qubit Hamiltonian with the listed parameters. 



on the Jc values derived from the IVCs of the readout SQUID and the measured junction 
areas. As shown in figure |3| there is very good agreement between the simulation result 
and the two-parameter fit for a = 0.6. However, significant deviations appear for 
a = 0.75. The reason is that there are not enough data points around 5$ = 0, where 
the readout of the qubit state by the dc SQUID fails. This leads to large uncertainties 
in A for the two-parameter fit. Therefore, small A values tend to have larger error bars. 
Nevertheless, figure |3] clearly demonstrates that the numerical simulation describes the 
experimental data very well. 

4^3. The fixed- gap gradiometric flux qubit 

We next discuss the properties of fixed-gap gradiometric fiux qubits to demonstrate 
the operation of the gradiometric qubit design shown in figures [4](a) and (b). The fiux 
qubit is biased close to its symmetry point by freezing in an odd number (2n + 1) of 
fiux quanta in the trapping loop during cool-down. This results in a phase difference 
of (2n + 1)7T between points A and B, equivalent to a fiux bias of (n + |)$o of the 
gradiometric fiux qubit at its symmetry point. To change the energy bias s after cool- 
down, a spatially inhomogeneous magnetic field is required, which is generated by the 
current sent through the 6:-fiux line. The qubit state is read out via the readout dc 
SQUID inductively coupled to the trapping loop of the qubit. The operation point of 
the readout dc SQUID can be optimized by applying a homogeneous magnetic field (e.g. 
by a solenoid) which does not affect the energy bias of the qubit due to its gradiometric 
design. 

Figure [41(c) shows typical spectroscopy data of a gradiometric fiux qubit with 
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ao = 0.65. Since we are measuring (jJc^{SI^) and not CL;q(5/), the only problem in 
evaluating this data is to determine the calibration factor 

where 61^ = — /|^"^ is the deviation of the current sent through the 6:-flux line 
from the value /|^"^ needed for biasing the qubit at the symmetry point. This is done 
by calculating uJq{Sf) by numerical simulations using Ec and a = A^/Aj as input 
parameters. The scaling factor tv is then obtained by re-scaling the measured uj^{61^) 
dependence to obtain optimum agreement with the simulation result. For the sample 
in figure |4](b), we obtain = 0.7mA~^, saying that a current of about 1mA results in 
Sf = 1. In general, the agreement between the experimental data and the simulation 
was found to be very good. The simulated values for the sample in figure [4](c) are 
A/27r = 5.1 GHz, /p = 420 nA and ao = 0.65. Again, we can make a consistency check 
by calculating the /p value according to ^ as discussed above. We obtain /p = 485 nA 
in good agreement with the value derived from the simulation. 

We note that we can also trap an even number 2n of fiux quanta in the trapping 
loop. In this case the phase difference between points A and B is 27m. This corresponds 
to a fiux bias of the gradiometric qubit by 2n$o/2 = n$o instead of (n+ |)$o for an odd 
number of trapped fiux quanta. That is, the qubit is biased far away from its symmetry 
point and no qubit transitions should be observable. This is in full agreement with 
the experimental observation. We finally note that I^r = {n — /ex)$o/(^g + ^k) is not 
allowed to exceed the critical current of the trapping loop. For our samples, this fact 
limits the number of trapped fiux quanta to a maximum value of 10 — 15. 

We also use the simple fixed-gap gradiometric qubit to check the quality of the 
gradiometer discussed in section 2.6 Figure [5](a) shows the switching current of the 
readout SQUID as a function of 6fi2 = fu — + |) [cf- (18)] recorded for a fixed 
microwave frequency of 19.33 GHz. The peaks and dips in the Iswi^fu) curves mark the 
^fi2 positions where the qubit transition frequency (jJq/27T = 19.33 GHz. On varying the 
number n of trapped fiux quanta, these positions shift due to the imperfect balance of 
the gradiometer. From the measured shift we derive (5grad,n = 943 ± 19. In figure [5](b), 
Sf 12 is plotted versus /ex generated by a homogeneous applied magnetic field. From 
the measured slope the quality factor (5grad,ex = 1076 ± 16 is determined. The total 
quality of the gradiometer is then Q 500, corresponding to a gradiometer imbalance 
of only 0.2%. This means, that the qubit operation point is shifted by about 2 m$o when 
we apply a homogeneous field generating one $o in the trapping loop. The measured 
quality factors are plausible. For example, the limited precision of the electron beam 
lithography process causes a finite precision 5A/Atr of the trapping loop area as well 
as 6S/S of the cross-sectional area and 6£/£ of the length of the superconducting lines. 
The measured quality factor corresponds to SA 0.2 /xm^, SS 50 nm^ or Si 60 nm. 
These values agree well with the values expected for the precision of the fabrication 
process. 

In figure [5|(b) we also plot the distance between the peak and dip positions in the 
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Figure 5. (a) Switching current I^w of the readout SQUID of a fixed-gap gradiometric 
flux qubit as a function of Sf 12 = /12 — + ^) recorded for a fixed microwave 
frequency of 19.33 GHz. The peak and dip positions mark those df values where 
co'q/27r = 19.33 GHz. (b) Frustration imbalance Sfi2 as a function of the frustration 
/ex generated by a homogeneous applied magnetic field. Also shown is the distance 
between the peak and dip positions in the Isw{^fi2) curves. 



/sw(5/i2) curves. This distance is about independent of /ex- This demonstrates that the 
qubit potential is not affected by the homogeneous background field. In total, our results 
show that the gradiometric fiux qubits can be fabricated in a controlled way and work 
as expected. The fact that the qubit operation point is not affected by a homogeneous 
background field allows us to integrate these qubits into large scale circuits where several 
qubits have to be operated and read out simultaneously without affecting each other. 



4^4' The tunable- gap gradiometric flux qubit 

In this subsection we discuss the results obtained with tunable-gap gradiometric fiux 
qubits as sketched in figure [l](d) . Besides a step- wise variation of a by freezing in an 
odd number of fiux quanta in the trapping loop we can make a continuous variation 
of a by an applied magnetic field generated either by the current /coii fed through an 
external solenoid or the current fed through the on-chip a-fiux line. We first discuss 
the experiments using a homogeneous magnetic field of a solenoid placed underneath 
the sample. The homogeneous magnetic field generates the frustrations /tr,net and fa,net 



of the trapping and a-loop, respectively, which are given by (11) and (12). 

Spectroscopy data of a tunable-gap gradiometric fiux qubit is shown in figure [6](c). 
The different a values were generated by the homogeneous magnetic field of the solenoid, 
whereas the fiux trapped during cool-down was constant at a single fiux quantum, i.e. 
n = 1. We can fit the data by a two-parameter fit yielding A and the slope duj^/d8I^ 
at large uj^ values. Here, 61^ = — /|^"^ is the deviation of the current sent through 
the 6:-fiux line from the value /|^"^ needed for biasing the qubit at the symmetry point. 
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loop 





Figure 6. (a) Circuit schematics of the tunable-gap gradiometric flux qubit with 
readout dc SQUID, e- and a-flux hue. The outer loop of the flux qubit (broken olive 
line) forms the trapping loop, the inner (grey line) the a-loop. (b) Scanning electron 
microscopy (SEM) image of the implemented circuit, (c) Transition frequency cjq/27r 
plotted versus 51^ = I^ — I^^^ for three different a values for a tunable-gap gradiometric 
flux qubit with ao = 1.10. Also shown is the result of a two-parameter flt. 



To derive the persistent current /p = {h/2^Q){d(jjJd5f) from this slope, we have to 
calibrate the horizontal axis. For this we need the calibration factor n = d6f /d6I^ [cf. 
(19)], which has already been discussed above. 

For the analysis of the A{a) dependence we need a second transfer function, relating 
the coil current /coii sent through the solenoid to the frustration fa^net of the a-loop. With 



(11) and (12) we obtain 

,net 



1 



(20) 



a/coil At, 1+ /3 a/coil ' 

With this transfer function and the expressions ^ and (13) for Ip and a, respectively, 



we obtain 



dujq 

dSf 



2$o/p 

h 



2$n/, 



O-'c 



h 



2«o 



cos VrC/coil + TT 



Ag /3 
Atrl+P' 



-n 



(21) 



Using the abbreviations rj = 2^oIcK/h and = (y4.o,/y4.tr)(/5/l + /3)(n/Q this simphiies 
to 



dSL 



r]Jl- [2ao I cos (7rC[/coii + /„ 



(22) 



We can use this expression to fit the measured duj^{Icoi\) / 051^ dependence using ry, 1^ 
and C as fitting parameters. 
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Figure 7. (a) Measured duo^/ddl^ values plotted versus the coil current IcoH through 
the solenoid producing the homogeneous magnetic field for a tunable-gap gradiometric 
fiux qubit. The solid line is a fit to the data using (22) yielding ( and In- The 
broken line shows the calculated a(/coii) dependence for these fitting parameters, (b) 
Minimal qubit transition frequency A/27r plotted versus a. The solid line is obtained 
by numerical simulations based on the full qubit Hamiltonian using the parameters 
Ej/h = 200 GHz and Ejh = 1.6 GHz. 



In figure [7|(a) the measured duj^/d5I^ values are plotted versus /coii together with a 
fit by (22). Evidently, the data points are clustered near specific /coii values. The reason 
is that the homogeneous magnetic field produced by /coii also changes the frustration of 
the readout SQUID and that the sensitivity of this SQUID is suflBcient only in a limited 
range of frustration. Figure [7|( a) shows that the expression (22) fits the experimental 
data well, yielding values for C, and With these fitting parameters we can calculate 
a = Q^o |cos (7rC[/coii + ^n])|- The resulting curve is also shown in figure [7|(a). We note, 
however, that in this case the fit parameters In and C cannot be used to directly 
determine /3 from the expression 1^ = (^a/Ar)(/3/l + /3)(n/C), because the value of 
In can be distorted by an additional background magnetic field. Therefore, we use 
only diflFerences A/^ to determine /3 [cf. (23)]. Knowing the Q^(/coii) dependence, we can 
adjust a to any desired value by adjusting /coii and then do spectroscopy at these values. 
Fitting the spectroscopy data [cf. figure [6](c)], we can derive the qubit gap A and plot 
it versus a. The result is shown in figure [7](b) together with the dependence obtained 
from numerical simulations based on the full Hamiltonian. The agreement between 
the experimental data and the numerical simulation is best for Ej/h = 200 GHz and 
Ec/h = 1.6 GHz, i.e. Ej/E^ = 125. We note that the Ej value agrees well with the one 
estimated independently from the measured junction areas and the Jc value measured 
for the junctions of the readout SQUID. This clearly shows the consistency of the data 
analysis and demonstrates the good control on the junction parameters fabricated on 
the same chip. Knowing the A(q{) and Q^(/coii) dependencies we can adjust the qubit gap 
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Figure 8. (a) Transition frequency a;q/27r plotted versus ^I^ = I^ — I^^^ for a tunable- 
gap gradiometric flux qubit with ao = 0.7. Also shown is the result of a two-parameter 
fit. (b) Measured duJc^/ddl^ values plotted versus the coil current /coil 

through the 

solenoid producing the homogeneous magnetic field for n = —3 trapped flux quanta. 



The solid line is a fit to the data by (22) yielding and 7^. The broken line shows the 
calculated a (/coil) dependence for these fitting parameters, (c) Measured values as in 
(b) for three different values of the number of trapped flux quanta, n = —3,-1, +3, 
fitted with consistent parameters. From the horizontal displacement of the different 
curves we obtain A/^/An = 0.43 mA. (d) Minimal qubit transition frequency A/27r 
plotted versus a for three different values of trapped flux quanta. The solid line 
is a fit of the data based on the full qubit Hamiltonian with the fitting parameters 
E:i/h = 800 GHz and Ejh = 1.0 GHz. 



in situ by /coii, while operating the qubit at the symmetry point with optimal coherence 
properties. This is a key prerequisite for many applications of flux qubits. 

For the sample of flgure [7| the qubit gap could be varied between values close 
to zero and about 5 GHz. For comparison, the data of a sample with larger ratio 
Ej/Ec = 800 is shown in flgure [sj The overall behaviour of this sample is very similar 
but the qubit gap can be tuned to values above 10 GHz. Moreover, we investigate the 
tunability of this sample for diflFerent amounts of trapped flux quanta n. In flgure [81(c) 
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we plot duj^/d5I^ versus the coil current /coii through the solenoid for three different 
values of the trapped flux ranging from n = — 3 to n = +3. Evidently, the general 
shape of the three curves is very similar as well as the obtained fltting parameters ( 
and 7]. The shift along the horizontal axis is expected from (22) and can now be used 
to calculate /3. Starting with the expression 1^ = {Aa/Atj.){/3/l + /3){n/()^ we only 
use differences A/^ = /^^^ — I^j. They correspond to differences An = rii — rij and 
result in A/^ = (Ac,/Ar)(/5/l + /3){An/(). For our sample, we flnd a mean value of 
A/^/An = 0.43 mA, flnally yielding 

An Ar ^ 



This value is in reasonable agreement with the one derived from the Lg and 
values which can be estimated from the qubit geometry, the cross-sectional area of 
the superconducting lines and the dirty limit expression of Lk. We note that the result 
from (23) is more precise because it is computed directly from the sample. All in all, 
our results show that the measured data agree well with the behaviour expected from 
theory. Moreover, the values of Ej and Ec obtained from fltting the data agree well with 
those obtained for junctions fabricated on the same chip. This demonstrates that the 
gap of gradiometric flux qubits can be reliably tuned over a wide range, making them 
attractive for a large number of applications. 

We flnally address the tuning of A by the on-chip a-flux line. Since the maximum 
current through this line is limited by its critical current and by heating effects in 
contacts, only small variations of the frustration of the a-loop are possible. Therefore, 
a constant applied magnetic fleld or a proper number of trapped flux quanta are used 
to pre-bias the qubit at a value o^b, where the slope of the A (a) dependence is steep. 
Then, 1^ is used to vary a around this value. In our experiments, a constant applied 
fleld is used to set o^b- Since the variation of the frustration of the a-loop is generated 
by la instead of /coii, we have to use the modifled calibration factor 

C^^^. (24) 

With this factor, we obtain 



cos ^arccos (ofb/ofo) + 7r(Ia^ . (25) 



051, 

We can use this expression to flt the measured du^/dSI^ versus 1^ dependence using 
C and Tj as fltting parameters. Based on these results, we can calculate a = 



cos ^arccos(ab) + 7t(Ic 



Knowing a, we can use the A values obtained from two- 
parameter flts of the spectroscopy data to get the A{a) dependence. Experimental data 
for the two samples of flgures [7] and [s] are shown in flgure |9| In flgures [9](b) and (d), 
we compare the experimental A{a) curves to numerical simulations based on the full 
qubit Hamiltonian with the same Ej and Ec values as obtained by tuning a with the coil 
current [cf. flgure[7|(b) and flgure[8](d)]. The very good agreement between measurement 
data and calculation demonstrates again the consistency of our data analysis. All in all. 
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Figure 9. (a) and (c) Measured duc^/dSI^ values plotted versus the current la 
through the a-flux line for two tunable-gap gradiometric flux qubits. The solid lines 
are fits to the data by (25) yielding the fitting parameters ( and rj. (b) and (d) Minimal 
qubit transition frequency A/27r plotted versus a of the qubits of (a) and (c). The 
solid lines are obtained from numerical simulations based on the full qubit Hamiltonian 
using the parameters Ej/h and E^/h as listed in the subfigures. The data in (a,b) and 
(c,d) are obtained for the two samples of figure [t] and figure [sj respectively, however 
with on-chip control of A via the a-flux line. 



our data clearly show that the qubit gap can be varied in a controlled way over a wide 
range by varying the frustration of the a-loop of the gradiometric flux qubits either by 
an external solenoid or an on-chip control line. 

5. Summary 

In summary, we have designed and fabricated gradiometric flux qubits with flxed and 
tunable gap. The characteristic parameters of the qubits have been derived from 
spectroscopy measurements. By trapping an odd number of flux quanta in the outer 
gradiometer loop during cool-down in a constant applied magnetic fleld, we were able 
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to pre-bias the gradiometric flux qubits at the symmetry point. We also performed a 
systematic analysis of the effect of the kinetic inductance of the narrow superconducting 
lines forming the qubit loop. The experimental results are in good agreement with the 
theoretically expected behaviour. The detailed analysis of the gradiometer imbalance 
showed that we can fabricate gradiometric qubits with an imbalance as small as 0.2%. 
This gradiometer quality is sufficient for most applications. 

Since the tunability of the qubit gap of persistent current qubits is a key issue, we 
have performed a systematic study on the tuning of the gap of gradiometric flux qubits 
by external control parameters. To this end, we have replaced one of the Josephson 
junctions in the qubit loop by a dc SQUID. This allowed us to tune the critical current 
of this junction and, in turn, the qubit gap in situ by a control flux threading the 
SQUID loop. The control flux was generated by three different methods: (i) an external 
solenoid, (ii) a persistent current frozen into the outer gradiometer loop, or (iii) a current 
sent through an on-chip control line. We have performed spectroscopic measurements, 
demonstrating a well-defined controllability of the qubit gap between values close to 
zero and more than 10 GHz. Our results clearly show that it is possible to vary the 
qubit gap over a wide range, as it is desired for tuning in and out of resonance with 
superconducting quantum circuits, while operating the qubit at its symmetry point 
with optimal dephasing properties. Due to the steep dependence of the qubit gap on 
the control ffux in some parameter regime, a very fast tuning of the qubit gap with 
small currents through on-chip control lines is feasible. 

We have compared the experimental data to model calculations based on the full 
qubit Hamiltonian. In general, very good agreement between experiment and model 
calculations is achieved. Fitting the data allowed us to determine the Josephson coupling 
and the charging energies of the qubit junctions. The derived values agree well with 
those measured for single junctions or SQUIDs fabricated on the same chip. By the 
controlled tunability of the ffux qubits a major drawback of this qubit type has been 
overcome. With their specific advantages such as their large anharmonicity and their 
potentially strong coupling to resonators, tunable-gap gradiometric fiux qubits are highly 
attractive for the implementation of quantum information circuits or the realization of 
fundamental quantum experiments. 
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